Optical Production of Stable Ultracold 88 Sr 2 Molecules 
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We have produced large samples of stable ultracold 88 Sr2 molecules in the electronic ground 
state in an optical lattice. The fast, all-optical method of molecule creation involves a near- 
intercombination-line photoassociation pulse followed by spontaneous emission with a near-unity 
Franck-Condon factor. The detection uses excitation to a weakly bound electronically excited vi- 
brational level corresponding to a very large dimer and yields a high-Q molecular vibronic resonance. 
This is the first of two steps needed to create deeply bound 88 Sr2 for frequency metrology and ul- 
tracold chemistry 

PACS numbers: 67.85.-d, 34.50.Rk, 37.10.Jk, 37.10.Pq 



The rapid progress in laser cooling has given rise to 
many new fields of research. One important example 
is the study of ultracold, dense clouds of molecules. 
The molecules can exhibit new physical phenomena near 
quantum degeneracy [H-Q. For example, long-range 
anisotropic interactions are expected between heteronu- 
clear polar molecules. Such molecules have also been ex- 
plored as a paradigm for quantum information and com- 
putation On the other hand, homonuclear molecu- 
lar dimers without a dipole moment present a metrolog- 
ical interest, for example in constraining the variation 
of the electron-proton mass ratio 0, ||| or complement- 
ing atomic clocks by serving as time standards in the 
terahertz regime Q . These molecules also provide an ex- 
cellent testing ground for many possible approaches to 
creating large ultracold samples, trapping them to al- 
low long interrogation times, and precisely controlling 
their quantum states. The four promising routes toward 
trapped neutral molecules [l[ are direct control of polar 
molecule dynamics via electric fields; buffer gas cooling 
of magnetic species; direct laser cooling of a suitable class 
of molecules [7j; and using magnetic or optical fields to 



combine laser cooled atoms into dimers [8j-[12j]. The lat- 
ter process typically results in molecules with relatively 
small binding energies but nonetheless has played a ma- 
jor role in the study of new phenomena. These dimers 
are also most promising for quantum control, since they 
are routinely produced at sub-/iK temperatures. 

In this Letter, we describe optical production of 88 Sr2 
in the electronic ground state in an optical lattice. In 
contrast to alkali-metal atoms, Sr atoms possess no elec- 
tronic spin and cannot be combined into molecules via 
the magnetic Feshbach resonance technique [13(. How- 
ever, Sr has the advantage of the intercombination (spin- 
forbidden) transition from the ground state, 1 Sq — 3 P\ 
(7 kHz linewidth [l4j, 689 nm wavelength), which al- 
lows Doppler cooling to < 1 /iK [l5[ and provides sev- 
eral features that enable efficient photoassociation (PA) 



logical applications. Furthermore, there is active inter- 
est in exploring molecules of alkaline-earth-metal (and 
the isoelectronic Yb) atoms in various combinations with 
each other and with alkali-met al atom s, such as Sr2, Yb2, 
LiYb, 
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into molecules [14j, llfjj. The unusually fast time scale 



(0.25 s) of ultracold molecule production is particularly 
important for establishing a short duty cycle in metro- 



RbYb, RbSr, and SrYb 
Our approach takes advantage of several interest- 
ing properties presented by the narrow intercombina- 
tion transition to achieve efficient transfer of atoms into 
molecules. Figure 03a ) and its caption explain the nota- 
tion used in this work. The excited molecular potential 
dissociating to the 1 Sq + 3 Pi atomic limit has a small 
attractive C3 coefficient, with a C§ coefficient similar to 
that of the ground state, resulting in relatively large wave 
function overlaps with the ground state van der Waals 
potential. This leads to several Franck-Condon factors 
(FCFs) between high- lying vibrational states that are un- 
usually close to unity. We exploit these large FCFs for 
one-photon PA that is followed by spontaneous decay to 
predominantly a single vibrational level. Another feature 
enabled by the narrow spectral line is the possibility to 
frequency resolve the least-bound electronically excited 
vibrational level with only a 440 hxkHz binding energy 
corresponding to the dimer size exceeding 500 ao [3] (h 
is the Planck constant and ao « 0.053 nm is the Bohr 
radius). This level has an optical length l opt (the cou- 
pling strength proportional to the PA laser intensity and 
the free-bound wave function overlap) that is estimated 
to be ~ 10 4 times larger than for the vibrational lev- 
els bound by hundreds of megahertz. We tune the atom 
recovery laser to the bound-bound transition resonant 
with this level, and a significant fraction of subsequent 
spontaneous decays leads to a reappearance of ground 
state atoms. Finally, this molecule creation method is 
the first step in a two-step sequence of producing Sr 2 
in the absolute ground quantum state. To this end, re- 
cent work [24], l25j predicts efficient two-photon coupling 
of the vibrational levels v = —2,-3 (Fig. [1]) to the low- 
est level v = via an intermediate v' ~ 12. We have 
made the molecules in v = —2 via PA into v' — —5 (3.5 
/ixGHz binding energy); furthermore, the calculated and 
measured FCF is also near unity for (v,v r ) = (—3,-6), 
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FIG. 1: (a) Long-range 88 Sr2 potential energies versus inter- 
nuclear separation. The two curves relevant to this work are 
the ground state potential and the excited state po- 

tential (l)0j that dissociate to the 1 So + 1 So and 1 So + 3 Pi 
(689 nm) atomic limits, respectively. Several vibrational en- 
ergy levels v,v', with the total molecular angular momenta 
(J, J') = (0, 1), are shown, along with their binding energies; 
the primed values refer to the excited potential. The negative 
v,v' values refer to counting from the dissociation limit, with 
— 1 corresponding to the highest-lying level. The natural line 
widths of the shown v' levels are ~ 20 kHz. The horizon- 
tal dashed line represents the thermal continuum for the fjK 
atoms. The vertical solid arrows indicate one-photon optical 
pathways used for molecule creation and atom recovery. The 
vertical dashed arrows indicate decay pathways completing 
the molecule creation and recovery; these rely on the unusu- 
ally large FCF and optical length, respectively, (b) The large 
value of the FCF f vv i — /(_2,-5) ~ 0.8 results from an excel- 
lent agreement of the outer turning points, as illustrated by 
the wave function plots for these vibrational levels. 



and our technique can be used for PA into v' — —6 (8.4 
hx GHz binding energy) by attaining a higher power out- 
put of the PA laser. The ability to make deeply bound 
Si'2 can enable a new class of highly precise molecular 
metrology tools [6|. 

Producing 88 Sr2 in the electronic ground state required 
performing precise two-photon PA spectroscopy. The 
search for v — — 2 and v — — 3 was guided by the most 
accurate ground and excited state potentials available 
from Refs. [H, [27| . Starting with 5 x 10 5 1 fiK 88 Sr 



FIG. 2: Autler-Townes doublets emerge when the Lbb fre- 
quency is fixed near a molecular resonance while Lfb is 
scanned, (a) Autler-Townes peak positions, along with their 
separation Q' fitted to Eq. JT]), for (v,v') = ( — 2,-5). 
The dashed line indicates the on-resonance Rabi frequency 
£1 at the indicated Lbb power, (b) The same as (a), for 
(v,v') — (—3, —4). The Lbb power is 60x larger, while Q is 
4.4x smaller, resulting in /(_ 3i _4)//(_2 l _5) = 0.9x 10~ 3 . The 
inset shows a representative atom loss curve with an Autler- 
Townes doublet corresponding to the circled data points. 



lattice formed at the 25 /mi waist of a retrorcflcctcd 914 
nm [28j laser beam (axial trap frequency 65 kHz; 95% 
of atoms in the zero-point trap level), we carried out the 
spectroscopy via intermediate electronically excited lev- 
els, with varying detunings from v' . The polarizations of 
the copropagating PA and detection lasers in this work 
are parallel to each other and to the small residual mag- 
netic field at the atom trapping site and perpendicular to 
the lattice polarization; the lasers are aligned along the 
tight-trapping axis. When the bound-bound laser (Lbb) 
is near resonance, and the frequency of the free-bound 
laser (Lfb, phase-locked to Lbb via a cavity-stabilized 
cooling laser) is scanned, Autler-Townes splitting of the 
PA resonance is observed via atom loss [Fig. d^b), in- 
set]. This splitting is equivalent to the generalized Rabi 
frequency of the bound-bound transition 



(1) 



atoms trapped in a 20 fiK deep one-dimensional optical 



where uj is the frequency of Lbb, wo corresponds to 
the bound-bound resonance, and Vl vv i is the Rabi fre- 
quency. Near the dissociation limit, £l vv ' is directly re- 
lated to the dimensionless FCF f vv > — \(v\v')\ 2 , since 
fl vv > — fl a \/2f vv >a, where a 2 — 1/3 is the rotational line 
strength for (J, J') = (0, 1) transitions, and the atomic 
Rabi frequency is fl a = T a ^J s/2 in terms of the atomic 
spontaneous decay rate T a and saturation parameter s. 

Different (v,v') pairs yield complementary FCFs; the 
pairs relevant to this work are (—2,-5) and (—3,-6). 
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FIG. 3: (a) Solid line: Fraction of the atom cloud detected af- 
ter a molecule producing PA pulse followed by an atom clear- 
ing pulse and an atom recovery pulse. The data are scaled to 
the calculation; the lowest bound on the peak atom-molecule 
conversion fraction is 3%. Dashed line: Calculated fraction of 
the atom cloud converted to v = —2 molecules after the PA 
pulse. The curve parameters were measured independently, 
with the exception of the molecule lifetime, which was fixed at 
3 ms to best match the data. The inset shows denser data in 
the first millisecond, (b) Fraction of the atom cloud detected 
after a 1 ms PA pulse followed by the atom clearing and 20 
fiW atom recovery pulses, as a function of the Lbb frequency. 
The sharp peak corresponds to the bound-bound transition 
from v = —2 to v' = — 1. Its fitted Lorentzian full width of 
32(4) kHz corresponds to a quality factor Q — 1.4 x 10 10 , the 
highest reported for a molecular vibronic transition, (c) The 
same as (b), at the higher 380 fiW atom recovery pulse power. 
In this high power regime, atoms are also recovered via the 
1 So + 3 Pi thermal continuum, which is only 0.4 MHz away 
from v = — 1. 



The corresponding FCFs were calculated to be 0.8(1), in- 
dicating a particularly strong coupling between the levels. 
Figure Eta) shows the Autler-Townes peak positions for 
(v, v') = (—2, —5) at various frequencies of Lbb- The dif- 



ference between the two curve branches is Q,', 



(-2,- 



and is 



-5) 

fitted to Eq. p. The FCF is found to be /(- 2 ,-5) = 1.1, 
which is consistent with 0.8(1) if the 27 /im PA laser beam 



waist measurement has an error of just 15%. In contrast, 
Fig. HJb) shows the Autler-Townes peak positions and fit 
for (v,v') = (—3,-4). With 60x more Lbb power, the 
splitting is 4.4 x smaller, leading to the experimentally 
determined ratio /(-3,_4)//(_2,-5) = 0.9 x 10 -3 . This 
agrees well with the calculated ratio ofl.lxl0~ 3 . While 
the (v, v') = (—3, —6) pair was not used in this work, our 
measurements have confirmed that its FCF is also near 
unity. Note that the fits such as shown in Fig. [2]yicld the 
binding energies for v = —2 and v = —3 levels (J = 0) 
as 1400.1(2)/ixMHz and 5110.6(2)/ixMHz, which agree 
with our calculations to 1.3% and 0.3%. 

The FCFs between weakly bound vibrational levels 
were computed from adiabatic single-channel potentials 
fitted to high resolution Fourier transform spectra of 
Refs. [H,[27|]. They were also calculated in Refs. [H[25| 
by describing the rovibrational dynamics in the excited 
electronic states by a coupled multichannel fully nona- 
diabatic Hamiltonian. The large calculated value of 
/(-2,-5) = 0.8(1) is surprising at first glance, since the 
wave function of v' — —5 shows 38 oscillations, while the 
wave function of v = —2 shows 61. However, according 
to the Franck-Condon principle [29], molecular transi- 
tions most readily occur between rovibrational levels of 
two different electronic states with nearly the same classi- 
cal turning points. In the adiabatic picture the classical 
turning points of v' = —5 and v = — 2 are very close, 
with their respective inner points at 7.1 and 7.5oo and 
outer points at 47.1 and 50.0ao- This near-coincidence 
explains the exceptionally large overlap between the two 
wave functions at large interatomic distances, as plotted 
in Fig. QJb). A similarly large overlap is predicted for 
v' = —6 and v = — 3, with the respective outer turning 
points at 40.5 and 40.1a [30| | . 

The 88 Sr2 molecule formation relies on one-photon PA 
into v' = —5 followed by spontaneous decay into v = —2. 
The PA induced atom density loss from the lattice is 
described by n a = —2Kn 2 a: where K is the two-body 
loss rate and we have neglected the finite lifetime of the 
trapped atoms (~ 10 s). Thus n a (t) = n /(l + 2Knot), 
where no is the initial atom density. Hence the time evo- 
lution of the molecule density is approximately given by 



0.3/w^< 



j/t c , where r c = (n 7 c ) is the 



lifetime parameter for the created molecules assuming 
similar rates for molecule-atom and molecule-molecule 
collisional losses, and the 0.3 factor results from the spon- 
taneous decay branching ratio from J' = 1 to J = 
[3lj |. (Note that with a simple repumping scheme out of 
J = 2, a tripled atom- molecule conversion efficiency is 
expected.) Therefore, the created molecule density as a 
fraction of hq can be expressed as 



d 
dt 



m 



n 
no 



(1 + atf 



m 



n 
no 



(2) 



For PA to v' = —5, the rates a = 2Kuq (and there- 
fore b = 0.3f vv 'Kno) were determined by measuring 
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atom loss as a function of the Lfb pulse length. At 
the Lfb power of 3.2 mW, we find a — 1680/s and 
b = 220/s. The loss rate c also has to be experimen- 
tally determined. If molecule-atom collisions dominate, 
it can be estimated as c ~ 2hnoRe/fJ- ~ 200/s [H} 
for our densities ^ ~ 2 x 10 12 /cm 3 , where /i is the 
molecule-atom reduced mass; the van der Waals length 
R 6 = 0.5(2 A iC* 6 /?i 2 ) 1 / 4 5.0 nm, where % = h/{2n) and 
Cq is the relevant attractive coefficient. This ~ 5 ms esti- 
mated collisional lifetime necessitates forming molecules 
on the millisecond time scale. While it requires relatively 
large intensities of Lfb , the undesirable photon scatter- 
ing by the atomic transition is suppressed due to the nar- 
row line. The numerical solution (2n m /no) to Eq. ([2]) is 
plotted in Fig. [3ja) for various PA pulse durations. As 
expected, for a sufficiently large PA rate a, the number of 
photoassociated atoms first rapidly grows and then grad- 
ually falls as the production slows relative to collisional 
loss. 

After the PA pulse is applied to the atoms, they are ex- 
posed to a low-intensity clearing pulse resonant with the 
strong 1 Sq — 1 Pi transition at 461 nm. The 0.2 ms pulse 
length is twice the minimum duration needed to remove 
the nonphotoassociated atoms to below the background 
level of the imaging setup; no systematic dependence on 
longer clearing pulse durations is observed. Then a 0.1 
ms atom recovery pulse is applied to the molecules cre- 
ated in v — —2. This pulse is resonant with (v, v') — 
(-2,-1); the calculated /(-2,-i) = 1-4 x 10~ 5 . As dis- 
cussed earlier, v' = — 1 has an exceptionally strong cou- 
pling to the free atom state: Z opt ~8x 10 5 ao/(W/cm 2 ) 
from calculations and ~ 4 x 10 5 ao/(W/cm 2 ) from pre- 
vious measurements [3]. Because of the small binding 
energy and large i op t, a significant fraction (30%) of the 
v' = —1 molecules is expected to spontaneously decay 
into ground state atoms that get recaptured in our opti- 
cal lattice. 

Figure [3(a) shows the result of the 88 Sr2 electronic 
ground state molecule creation followed by conversion of 
the molecules back to atoms. An overall scaling factor 
has been applied to the data. The expected curve is ob- 
tained from Eq. ([2]), with the molecule lifetime chosen 
to best match the data. This results in c ~ 380/s, or a 
3 ms lifetime of the v — —2 molecules. The calculated 
curve predicts that at the given Lfb power, 12% of the 
atoms should be converted to molecules at optimal PA 
pulse lengths of 1 ms. The actual fraction of recovered 
atoms is 1%, which corresponds to a conversion of 3% of 
the atoms into molecules assuming a 30% efficiency of the 
recovery pulse. Given this uncertainty, we conclude that 
we create 2(1) x 10 4 molecules with the 1 ms PA pulse. 
The lower-than-expected number of recovered atoms is 
likely due to a reduced efficiency of the atom recovery 
step rather than of molecule formation, since recapture 
depends on the lattice dynamics of atoms with a wide 
kinetic energy spectrum. Note that no significant heat- 



ing of the recovered atoms is observed at the end of the 
process. Detecting v = —2 molecules produced in the 
J = 2 state could yield more information on the overall 
efficiency. Furthermore, preliminary work indicates that 
purified 88 Sr2 samples can have very long lifetimes on the 
10 ms and 1 s scales in one- and three-dimensional lat- 
tices, respectively, limited by bimolecular collisions. De- 
tailed studies of molecule losses at these extended times 
in the regime of few molecules per lattice site could help 
ascertain the efficiency factors. 

Atom recovery via the least-bound excited molecular 
level is confirmed by choosing the optimal 1 ms molecule 
producing pulse duration and measuring the recovered 
atom fraction as a function of the Lbb frequency, as 
shown in Fig. [3][b). The narrow peak corresponds to 
the bound-bound transition from v = — 2 to v' = — 1. Its 
fitted Lorentzian full width of 32(4) kHz, approximately 
twice the expected natural width, corresponds to a high 
quality factor Q — 1.4 x 10 10 of the 88 Sr2 molecular vi- 
bronic transition. At stronger saturation, the resonance 
broadens as in Fig. (He), and the atoms can also be re- 
covered via the electronically excited thermal continuum, 
which is 0.4 MHz blue-detuned from the molecular line. 

In conclusion, we have produced stable ultracold 
lattice-confined 88 Sr2 molecules in the electronic ground 
state that are bound by 0.05 cm -1 . Our technique is 
the first step of a two-step projected sequence for achiev- 
ing large samples of ultracold Sr2 in the absolute ground 
quantum state. This has immediate applications in pre- 
cise time and frequency metrology, studies of fundamen- 
tal constant variations, and investigations of chemical re- 
actions at ultralow kinetic energies and highly nonther- 
mal internal state distributions. The all-optical molecule 
creation scheme relies on a large wave function over- 
lap between a pair of electronically excited and ground 
rovibrational levels that is characteristic of the alkaline- 
earth-metal-atom system. The entire experimental tra- 
jectory ending in the molecule sample takes only 0.25 
s. The atom recovery scheme, proving the creation of 
molecules, relies on one-photon excitation to the least- 
bound vibrational level near the intercombination line 
dissociation limit. This level can be resolved due to the 
narrow optical line width and has an unusually large 
coupling to the thermal continuum of the ground elec- 
tronic state. To achieve the goals of this work, we have 
performed one- and two-photon spectroscopic studies of 
88 Sr2. Note that related results are reported for 84 Sr in 
a Mott-insulator state [33l ]. 
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